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ABSTRACT 

Quasi-periodic oscillations and X-ray spectroscopy are powerful probes of 
black hole masses and accretion disks, and here we apply these diagnostics to 
an ultraluminous X-ray source (ULX) in the spiral galaxy NGC628 (M74). 
This object was observed four times over two years with the Chandra X- 
ray Observatory and XMM-Newton, with three long observations showing 
dramatic variability, distinguished by a series of outbursts with a quasi- 
period (QPO) of 4,000-7,000 seconds. This is unique behavior among both 
ULXs and Galactic X-ray binaries due to the combination of its burst-like 
peaks and deep troughs, its long quasi-periods, its high variation amplitudes 
of > 90%, and its substantial variability between observations. The X-ray 
spectra is fitted by an absorbed accretion disk plus a power-law component, 
suggesting the ULX was in a spectral state analogous to the Low Hard state 
or the Very High state of Galactic black hole X-ray binaries. A black hole 
mass of ~ 2-20 x 10 3 M Q is estimated from the fb~M, scaling relation found 
in the Galactic X-ray binaries and active galactic nuclei. 

Subject headings: Galaxy: individual(NGC628) — X-rays: binaries 

1. Introduction 

Ultra-luminous X-ray sources (ULXs), found in many nearby galaxies by EIN- 
STEIN, ROSAT, and recently by Chandra and XMM, are non-nuclear X-ray point 
sources with X-ray luminosities in the range of 10 39 - 10 41 erg/sec. If ULXs emit at 
the Eddington limit, such luminosities imply black holes (BHs) of masses 10-1000 M , 
or 100-10,000 M if the sources emit at 10% level of the Eddington luminosity. Such 
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intermediate mass black holes (IMBHs) can naturally explain the observed high lumi- 
nosities and, if they exist, bridge the gap between stellar mass BHs (~ 10 1 M Q ) and 
supermassive BHs (10 6 — 10 9 M Q ) in the center of galaxies. However, such IMBHs are 
not predicted to be the products of ordinary stellar evolution models, and it is still in 
debate whether they can form in dense stellar fields via runaway stellar collisions with 
seed BHs of a few hundred M (e.g., Portegies Zwart 2002). Alternatively, these sources 
can be stellar mass black holes, for which the formation is well established in theory 
and observation. Special mechanisms are required for these stellar mass black holes to 
emit at super-Eddington luminosities. Such a mechanism could be extreme relativis- 
ts beaming (Georganopoulos et al. 2002), mild geometric beaming (e.g., King et al. 
2001; Abramowics et al. 1978), or the photon-bubble instability in a radiation pressure 
dominated accretion disk, which leads to truly super-Eddington luminosities (Begelman 
2002). 

X-ray timing and spectral analyses have been widely used in the study of bright 
Galactic X-ray binaries. Combined with X-ray spectroscopy, rapid variability provides 
direct information on the compact object and the accretion disk, the site of X-ray emis- 
sion. Most Galactic X-ray binaries show variability characterized as periodic variations, 
quasi-periodic oscillations, and low frequency noise on times scales from 10~ 4 seconds 
to years. They emit in well-defined spectral states including the Very High (VH) state, 
the high/soft (HS) state, the low/hard (LH) state, and the quiescent state, (e.g., Lewin, 
Van Paradijs & Van den Heuvel 1995, and references therein). In comparison, ULXs are 
less variable on short times scales. Swartz et al. (2004) found in a recent Chandra sur- 
vey that only 5%-15% of the ULXs showed significant variability during their Chandra 
exposures. In another Chandra survey, Roberts et al. (2003) detected only short term 
flux variability in one of the two epochs for two out of the eight ULXs in their survey. 

Periodic variations and quasi-periodic oscillations are powerful probes of the orbit 
motion and the BH - accretion disk interaction, but have been detected only for a few 
ULXs so far. A ULX in IC342 showed small amplitude (5%) periodic variations on a 
time scale of 30-40 hours in a 155 hour ASCA observation (Sugiho et al. 2001). A 
ULX in the Circinus galaxy showed X-ray eclipses and a period of 7.5 hours, though 
it is possibly a foreground interloper in our galaxy given its proximity to the Galactic 
plane and the similarity to a long period AM Her system (Bauer et al. 2001). A ULX 
in M51 showed ~70% variations in a smooth sinusoidal light curve with a period of 
7620±500 seconds in a 15 ks Chandra observation (Liu et al. 2002), while in another 
27 ks observation, the luminosity dropped by a factor of 50 and there were too few 
counts to verify this periodicity. From XMM and RXTE observations of a ULX in M82, 
Strohmayer & Mushotzky (2003) discovered a narrow QPO of 54.4±0.9 mHz (~ 18 
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seconds) with a r.m.s. amplitude of 8.4%, which suggests a geometrically thin accretion 
flow and argues against substantial beaming. Soria et al. (2004) found that the light 
curve of a ULX in NGC5408 showed flares of amplitude of ~30% at a quasi-period of a 
few hundred seconds, which they attributed to variability in the power-law component 
of the emission. A ULX in the colliding galaxy NGC7714 showed medium amplitude 
(~30%) periodic variations on a time scale of 2 hours in a 15 ks XMM observation 
when it was in a high/hard state with L x = 6.6 x 10 40 erg/sec, while such a period 
was not detected in an earlier XMM observation when it was in a low/soft state with 
L x = 4.4 x 10 40 erg/sec (Soria & Motch 2004). 

A ULX in NGC628 (M74), CXOU J013651. 1 + 154547 was found to be extremely 
variable with unique quasi-periodic variation patterns in recent Chandra and XMM 
observations (Krauss et al. 2003). This source is located in a bubble nebula on a spiral 
arm of NGC628, a face-on Sc spiral galaxy with a galactic latitude of -45° at a distance of 
9.7 Mpc; it is most likely in NGC628 given the association with the nebula and absence 
of bright foreground objects or known AGNs/QSOs in the vicinity. In this paper, we 
analyze two Chandra observations and two XMM observations available in the archive, 
and report its spectral and timing properties in section 2. A discussion on the QPOs 
and the BH mass is presented in section 3. 



2. Data Analysis and Results 

There are four X-ray observations in the archive for the ULX, including two Chan- 
dra ACIS observations, on 2001 June 19 (Observation ID 2057, 46.3 kiloseconds) and 
2001 October 19 (ObsID 2058, 46.3 kiloseconds), and two XMM EPIC observations, 
on 2002 February 2 (ObsID 0154350101, 36.9 kiloseconds) and 2003 January 7 (ObsID 
0154350201, 24.9 kiloseconds). To extract spectra and event lists for the ULX, we used 
CIAO 2.3 and CALDB 2.20 for the Chandra data, and SAS 6.0.0 for the XMM data. 

For timing analysis, the event lists for the ULX in two Chandra observations were 
extracted from the 3a source regions enclosing 95% of the source photons as reported 
by the detection task WAVDETECT. To estimate the background photon counts in the 
source region, we extracted the events in a nearby region similar in shape but offset 
from any point sources. It was found that the background counts are less than 3% of 
the total counts. In the two XMM EPIC observations, the ULX fell into the chip gap 
on the pn, so only the two MOS observations are used for analysis. The events were 
extracted from a source circle with a radius of 512 pixels (i.e., 25'.'6) enclosing ~90% of 
the source photons, within which no other point sources are detected in the two Chandra 
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observations. We took the photons within a nearby circle offset from all point sources as 
an estimate of the background photons in the source circle; for the first (second) XMM 
observation, the background counts are estimated to be 34% (30%) of the total counts 
in the source circle. The light curves for the four observations were binned and plotted 
in Figure 1. 

Some distinct features of the light curves are immediately obvious in Figure 1. 
The burst-like peaks and deep troughs alternate quasi-periodically in the first three 
observations; the signal is too weak to see such features in the 4th observation. The 
changes in fluxes are violent on time scales of hours. In ObsID 2057, there was a period 
of 5 kiloseconds without any photons detected, preceded by a decrease in count rates, and 
followed by an increase that leads to a large burst. The average count rate is 7.5 (20.6) 
ACIS count/ks for the 1st (2nd) Chandra observation, and 27.1 (15.0) MOS count/ks 
for the 3rd (4th) XMM observation. The standard deviations from the average count 
rates are 11.2/19.8/28.3/12.8 count/ks respectively, indicating variation amplitudes of 
90% or higher. 

To study the quasi-periodicity for the ULX, the power-density spectra (PDS) were 
calculated with XRONOS 5.19 (Figure 2) for the first three observations. The PDS for 
ObsID 2057 shows a peak indicative of a QPO around / ~2.5xl0 -4 Hz, which appears 
as a hump in an otherwise declining power-law (f~ a with a~l). This PDS resembles 
that of a Galactic BH binary in its Very High state or steep power-law (SPL) state as 
described by McClintock & Remillard (2003). However, due to the large errors and low 
resolution at low frequencies, this PDS also could be viewed as a flat-top spectrum at 
low frequencies that breaks into a declining power-law with a~1.5 around /& ~2.5xl0~ 4 
Hz, thus resembling that of a BH binary in its LH state or the hard state (McClintock & 
Remillard, 2003). The PDS for ObsID 2058 shows a significant broad peak that is more 
clearly interpreted as a QPO within the frequency range of 1-4 xl0~ 4 Hz. The PDS for 
the first XMM observation is similar to that of a BH binary in its LH state, with a break 
frequency (i.e., QPO frequency) around / ~1.5xl0~ 4 Hz and a declining power-law with 
a~1.2. No significant features are found in the PDS for the second XMM observation, 
partly due to its low count rate and short exposure. 

A wavelet technique (e.g., Torrence & Compo 1998; Hughes et al. 1998) is used to 
search for the quasi-periods for the QPOs and has certain advantages over the Fourier ap- 
proach. The conventional Fourier transform may not work well to reveal quasi-periodic 
signals at low frequencies, because the power associated with the data's window can 
appear at low frequencies, and non-periodic outbursts will spread power across the spec- 
trum. In comparison, the wavelet analysis studies the behavior of the signal on different 
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time scales as a function of time, and preserves the temporal locality of the signal. The 
global power spectra calculated from summing up the power spectra at all times are 
shown in the period space to facilitate the search (Figure 3). To assess the significance 
of a quasi-period in the global power spectrum of a signal, its power is compared to the 
distribution of the power spectra of both white and red noise random processes at that 
time scale (Torrence & Compo 1998). Three quasi-periods at ~ 4000 seconds, ~ 6000 
seconds and ~ 7000 seconds are found in the first three observations above the 99% 
significance level, i.e., the powers at the quasi-periods are higher than 99% of the red 
noise random processes. The quasi-period at ~ 7000 seconds (i.e., two hour) was present 
in all three observations. 

The spectra for the first three observations were fitted with three models, the power- 
law model (i.e., wabs +powerlaw) , the multi-color disk model (i.e., wabs -l-diskbb), and 
the disk plus power-law model (i.e., wabs -l-diskbb +powerlaw) using XSPEC 11.2.0; no 
spectrum was fit for the 4th observation due to problems with the Observational Data 
File (ODF). The spectra were binned to a minimum of 15 photons per bin, and fitted 
in the 0.3-8 keV band (Table 1). The best model for both Chandra observations is the 
disk plus power-law model. Despite the flux change by a factor of ~3, the two spectra 
are strikingly similar, with similar model parameters (i.e., the same disk temperature 
of ~0.2 keV, and the same power-law photon index ~1.5), and similar fractional fluxes 
from the disk component (~l/4) and the power-law component (~3/4). The spectral 
decomposition suggests the ULX might be in a spectral state analogous to the Low 
Hard state of Galactic BH binaries, or the Very High state given the presence of low 
frequency QPOs and the power-law photon index up to ~2.2 allowed by the fitting 
errors. The power-law models are also acceptable for the two Chandra observations, 
with similar best fit column density and photon index (~2.2). Three emission lines 
in the two Chandra observations were detected at ~ la significance levels, and the 
details will be studied in a future work. For the combined MOS spectrum in the third 
observation, the disk plus power-law model is unphysical since the power-law component 
increases with photon energy (r ~ —2.5); the power-law model is acceptable, and has 
the same model parameters as in two Chandra observations within the errors. 

For ObsID 2058, two separate spectra for the high and low flux states, 2058-H and 
2058-L, are extracted from the time intervals above and below 20 count/ks, a sensible 
but non-unique division point. The two spectra sample the photons in the bursts and 
troughs respectively, and are fitted with the three models as the above. Both the power- 
law model and the disk plus power-law model are acceptable for these two spectra; in 
both models the parameters (column density, photon index, and disk temperature) for 
the high flux spectrum are the same as those for the low flux spectrum. Such a similarity 
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between these two spectra is expected only if the QPOs are present in both the disk 
component and the power-law component of the X-ray emission. In accordance to this 
implication, the QPOs with similar variation amplitudes are seen in the light curves 
of low-energy photons (< 1 keV, dominated by the disk component) and high-energy 
photons (> 1 keV, dominated by the power-law component). 

3. Discussion 

QPOs are seen in the X-ray observations for CXOU J013651. 1+154547, one of the 
four ULXs with possible QPOs. The quasi-periodicity is directly visible in the light 
curves, and it is present in observations spanning 2 years with flux changes by a factor 
of three. This is verified by the timing analysis with two different techniques, which have 
revealed significant QPOs as peaks and/or breaks in the PDS. Quasi-periods of ~ 4-7 
kiloseconds are found from the observations, most evident with a wavelet decomposition 
technique. In both Chandra observations, the X-ray spectra are best fitted with a 
multi-color disk plus a power-law component with similar parameters, in which the disk 
components (T in ~ 0.2 keV) dominate the soft band and contribute a quarter of the 
total flux, while the power-law components dominate the hard band and contribute 
three quarters of the total flux. This spectral shape, combined with the presence of low 
frequency QPOs, suggests the ULX was in a spectral state analogous to the Low Hard 
state or the Very High state of Galactic BH binaries. 

Compared to QPOs in other ULXs and Galactic BH binaries, the QPO for CXOU 
J013651. 1 + 154547 is unique for the combination of its burst-like peaks and deep troughs, 
its long quasi-periods, and its extreme variability. The light curves of CXOU J013651. 1+154547 
exhibit substantial variability from quasi-period to quasi-period. Within a quasi-period 
of about two hours, the troughs are usually 10% - 20% of the average flux or lower, and 
the peaks are up to 5 times the average flux, reaching ~ 5x 10 39 erg/sec. In comparison, 
QPOs for other ULXs and Galactic X-ray binaries usually have smooth (sometimes si- 
nusoidal) quasi-periodic patterns and r.m.s. amplitudes of a few percent reaching up to 
30% only in rare occasions. The QPOs are present in both the accretion disk component 
and the power-law component. As a comparison, the QPOs for an ULX in NGC 5408 
with a quasi-period of a few hundred seconds are significant only in the hard power-law 
component (Soria et al. 2004). 

The BH mass can be estimated with the fb~M, scaling relation if the QPOs originate 
from the accretion disk as for Galactic X-ray binaries. The scaling relation between the 
QPO frequency and the BH mass has been found to exist in Galactic BH X-ray binaries 
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and active galactic nuclei, albeit with large uncertainties in the calibration (e.g., Belloni 
& Hasinger 1990; Markowitz et al. 2003). The PDS for CXOU J013651. 1+154547 is 
comparable to that of Cygnus X-l (~ 10M Q ) which at low frequencies exhibits a flat 
top that breaks at fbi of 0.04-0.4 Hz into a declining power-law with slope ~ —1 that 
breaks again at /&2 of a few Hz into a steeper power- law with slope ~ —2. If the 
observed QPO frequency ~ 2 x 10" 4 Hz corresponds to fbi as suggested by the PDS 
shape, the BH mass is estimated as ~2- 20xl0 3 M Q . The BH mass would be ~ 1O 5 M 
if the QPO frequency corresponds to /j, 2 , as is possible owing to the uncertainties in the 
PDS of these short observations. The above mass estimates are not conclusive, since 
the low frequency QPOs occur at different locations of the accretion disk for different 
accretion rates, thus their frequencies are not determined only by the BH mass. This 
mass estimate is different from the mass estimated from its luminosities (< 100M o ) if 
the system emits isotropically at above 10% of the Eddington limit as indicated by the 
emission state possibly analogous to the Very High state in Galactic BH binaries. 

More convincing and stringent mass estimates can be made with high frequency 
QPOs, especially the 3:2 twin-peak QPOs that appear as kHz QPOs in Galactic X-ray 
binaries. This 3:2 twin-peak QPO is thought to occur at a specific resonance radius in 
the inner region of the accretion disk that is fixed in terms of the gravitational radius 
of the BH. Thus its frequency is expected to scale with the BH mass. Abramowics et 
al. (2004) advanced a stringent scaling relation / = 2.8/M. kHz for such QPOs, which 
is verified in three microquasars and the central BH of our Galaxy, Sgr A*. If the 3:2 
twin-peak QPO is to occur in NGC628-ULX, we expect a frequency of ~ 0.3 Hz given 
the BH mass of ~ 10 4 M Q as inferred from the observed low frequency QPOs. However, 
to verify such a QPO would require a lengthy observation and the object would need to 
be brighter than usual. 
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gratefully acknowledge support for this work from NASA under grants HST-GO-09073. 
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Fig. 1.— The light curves for CXOU J013651. 1+154547 in four observations. Over- 
plotted for comparison are sinusoidal curves with the periods of 4000/4000/7000/6000 
seconds respectively. 
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Fig. 2. — The power-density spectra calculated with the Fourier transform for CXOU 
J013651. 1 + 154547 in the first three observations. A broad peak around l-4xl0~ 4 Hz is 
clearly present in the second Chandra observation that can be interpreted as a QPO. 
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Fig. 3. — The global power spectra calculated with a wavelet technique for CXOU 
J013651. 1+154547. The dotted (dashed) line indicates the significance level at which the 
power of a (quasi-)periodic signal is above 99% of the random realizations of a red (white) 
noise process. Above the 99% significance level of a red noise process are quasi-periods 
at ~4000 seconds and ~7000 seconds in the first observation, quasi-periods at ~4000 
seconds, ~6000 seconds and ~7000 seconds blended together in the second observation, 
and quasi-periods at ~7000 seconds and ~ 12000 seconds in the third observation. The 
7 kilosecond quasi-period can be seen in all three observations. 



- 12 - 



Table I. Spectral fits for CXOU J013651. 1+154547 



spectrum 


model 


N H 


r 


kT in 


xl/doi 


Prob 


L X (0.3-8) 


id^fcbb (0.3-8) 


L pl (0.3-8) 






(10 20 cm,- 2 ) 




(Kcv) 




% 


(10 38 erg/s) 


(10 38 erg/s) 


(10 38 erg/s) 




wabs(pl) 


6.4±4.6 


2.4±0.3 




0.87/50 


73.4 


3.8 




3.8 


acis2057 


wabs(diskbb) 


0.0±100 




0.46±0.03 


1.12/50 


26.5 


2.8 


2.8 






wabs(pl+diskbb) 


8.8±10.4 


1.5±0.7 


0.21±0.08 


0.81/48 


82.7 


4.5 


1.3 


3.2 




wabs(pl) 


6.2±2.2 


2.2±0.1 




1.10/109 


22.3 


11.6 




11.6 


acis2058 


wabs(diskbb) 


0±100 




0.64±0.03 


1.53/109 





8.4 


8.4 






wabs(pl+diskbb) 


7.2±4.6 


1.4±0.4 


0.24±0.05 


1.04/107 


36.4 


13.4 


3.3 


10.1 




wabs(pl) 


15.8±3.2 


2.1±0.2 




0.92/86 


67.6 


21.4 




21.4 


2058-H 


wabs(diskbb) 


2.9±1.7 




0.93±0.09 


1.11/86 


22.7 


17.8 


17.8 






wabs(diskbb+pl) 


28.7±13.8 


1.8±0.3 


0.16±0.03 


0.88/84 


78 


23.4 


3.3 


20.1 




wabs(pl) 


9.6±4.4 


2.3±0.3 




1.21/43 


16.7 


3.8 




3.8 


2058-L 


wabs(diskbb) 


0±5.3 




0.74±0.10 


1.44/43 


3.0 


3.2 


3.2 






wabs(pl+diskbb) 


9.1±9.1 


1.7±0.9 


0.25±0.13 


1.22/41 


15.1 


4.2 


0.9 


3.2 




wabs(pl) 


14.3±8.9 


2.1±0.4 




1.01/35 


45.7 


4.6 




4.6 


0154350101 


wabs(diskbb) 


Oil 1.3 




0.91±0.17 


1.00/35 


46.7 


3.9 


3.9 






wabs(diskbb+pl) 


0±100 


-2.5±100 


0.86±0.39 


1.06/33 


38.1 


5.3 


3.8 


1.5 



